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We report on the design and experimental characterization of a surface-electrode multi¬ 
pole ion trap. Individual microscopic sugar particles are confined in the trap. The trajectories 
of driven particle motion are compared with a theoretical model, both to verify qualitative 
predictions of the model, and to measure the charge-to-mass ratio of the confined particle. 
The generation of harmonics of the driving frequency is observed as a key signature of the 
nonlinear nature of the trap. We remark on possible applications of our traps, including to 
mass spectrometry. 


I. INTRODUCTION 

Recent years have seen a wealth of innovation in the design of radiofrequency (RF) ion traps 
MM- To some extent, this has been driven by the utility of microfabrieated ion traps for quantum 
information processing (QIP). However, there is still aetive interest in the use of RF ion traps in 
atomic and molecular spectroscopy IfTTUTSlI . mass spectrometry [|T^ . nonlinear dynamics [fTTI - 
[T9]| . and other areas. Although most work geared toward QIP has foeused on quadrupole ion 
traps, whieh have an electric potential and time-averaged pseudopotential that are quadratic near 
the trap eenter, there are also groups pursuing new designs and applieations of multipole ion traps 
lHH l20] - l2^ . Our working definition of a multipole trap is that it is a trap that features a potential 
landscape that is of a higher order (going as or x®, for example, where x is the position of a 
trapped ion with respeet to the trap eenter). The reeent growth of interest in nonlinear dynamics 
using ions, in particular, has involved the use of quadrupole ion traps, but within whieh the ion is 
exeited to a high enough energy that it is acted on by higher-order terms in the trapping potential 
mm- For such experiments, a multipole trap would be advantageous. 

We have reeently proposed a trap design for a surface-eleetrode multipole ion trap (SEMT) 
t2M . There are a number of reasons why this design might be of interest to researehers. Firstly, 



the trapping potential is highly nonlinear near the trap eenter, whieh eliminates the need to impart 
as mueh energy to the ion to see nonlinear behavior, if that is one’s goal. Seeondly, no DC 
voltages are required for ion eonfinement, as is the ease with linear multipole traps. This results 
in less mieromotion than one would otherwise eneounter, whieh is benefieial in preeision laser 
speetroscopy and other applications. (Mieromotion is driven motion at the RF frequency which 
is present whenever the ion is located anywhere the RF driving field does not vanish. In linear 
multipole traps, the DC voltages on the endcaps result in electric fields that push ions further 
from the RF null than they would otherwise be [|22l .f Thirdly, the designs set forth in Ref. [[24ll 
are already being used, in a modified form, in the design of ion traps which will be at the heart of 
exotic experiments on time crystals [l25ll26l . Fourthly, ideas for new applications continue to arise, 
including the possibility of using SEMT’s as the basis of a new method of mass spectrometry. This 
will be discussed at the conclusion of this article. 

In this article, we describe the design and experimental characterization of a SEMT. Our key 
result is the demonstration of frequency upconversion in the trajectory of an ion which is close to 
the trap center, which is a hallmark of nonlinear behavior. In Sec. II, we describe the mathematical 
model of the trap. In Sec. Ill, we present our experimental setup and main results. Sec. IV contains 
a discussion of our results, and some reflections on possible applications of this kind of trap. 


II. TRAP DESIGN AND MODELING 

The design of our SEMT is presented in Eig. The trap was designed using the EAGEE 
PCB layout software and manufactured by Sunstone Circuits. It employs a silver finish to avoid 
corrosion and to make the electrodes equipotential surfaces to the greatest extent possible. The 
electrode radii chosen are only one of an infinite number of possibilities that will generate a 
multipole trap. These particular values were chosen because they roughly maximize the expected 
trap depth for a given RE voltage and frequency, and for an ion charge-to-mass ratio of 1 mC/kg 
(an order of magnitude estimate for dust particles). A surface-electrode multipole trap is created 
by applying alternating voltage to REl and RES while grounding RE2. 

The time-averaged effective potential which governs the motion of the ions, called the pseu¬ 
dopotential and denoted \l'(r, z), is given by 
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FIG. 1: Image of the trap as rendered by the EAGLE program. The overall width of the trap is 40.2 mm. 
The electrodes which may be set to a static voltage are labeled as “DCl” through “DC4.” The electrodes to 
which an oscillating voltage may be applied are labeled “RE2” and “RES.” The center electrode is “REl,” 
although it is not labeled. The nominal radii of the RE electrodes are 1.02 mm, 3.48 mm, and 8.62 mm. 
There is 0.15 mm gap between these electrodes. 


where Q is the eharge of the ion, m is the mass of the ion, G is the (angular) frequeney at whieh 
the trap is driven, and V<f>(r, z) is the RF field amplitude. We define the trap axis as the line 
perpendieular to the plane of the trap and intersecting the center of the trap electrodes. The z 
coordinate measures the displacement above the electrodes, measured along the trap axis, and 
the coordinate r measures the displacement perpendicular to the trap axis in a plane parallel to 
the trap electrodes. The calculation of <h(r, z) for many trap geometries, including ring traps, is 
a problem in electrostatics to which several authors have recently made important contributions 
ll27l - [30ll . In the present work, we use the methods outlined by Schmied ll30l . The purpose of using 
these methods is to take into account the effect of gaps between the electrodes, which were set to 
be the minimum allowed by the manufacturer (Sunstone Circuits). However, we do note that the 
true potential is indeed very close (within 1% inside the region of interest) to that obtained from 
the gapless approximation. In the gapless approximation, the electric potential is given by 
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$(^,r)= / Jo{kr)e-^^Ao{k)dk 


where Jq is the Bessel funetion of zeroth order and Ao^k) = The Ai are given by 


Ai{k) = Vi (ri+iJi(kri+i) - riJi(kri )), (3) 

where Vi is the amplitude of the RF voltage applied to the eleetrode, J\ is the Bessel funetion 
of first order, and the r* are the radii deseribing the trap geometry (as given in Fig. [^. There is 
a simple expression for the potential due to a given eleetrode earrying a voltage V along the trap 
axis (at r = 0), given by 


4(j) 
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where r 2 and ri are the inner and outer radii, respectively, for that particular electrode JHl . The 
full potential may be obtained by adding all potentials due to each RF electrode. In the work that 
follows, we focus on motion in the z direction; therefore, we will rely heavily on this expression. 

We now present the results of our modeling of the trap. As in the experiments described 
herein, we assume that the RF voltage applied to electrodes RFl and RF3 has an amplitude of 
660 V at a frequency of 297 Hz. All other electrodes are assumed to be at RF ground. We use 
the experimentally-obtained value of R = Q/m Ri —2.5 x 10“^ C/kg for the charge-to-mass ratio 
of a particular ion. A two-dimensional cross-section of \['(r, z) is shown in Fig. In Fig. we 
plot \l'(r = 0, z)/m and a fit thereof, in order to show the importance of the various higher-order 
terms. 

In experiments with atomic ions, the gravitational force on the particle can normally be ne¬ 
glected in calculating the pseudopotential. However, in the case of trapping micron-scale objects, 
the force of gravity is not negligible. This is especially so given the flat landscape of the RF 
potential near the trap center. In order to compensate the trap (which means to null the micro¬ 
motion of the ion), a DC voltage Vdc is applied to the electrode RF2. This DC voltage pushes 
the ion toward the point at which the RF electric field vanishes. For our present work, a value of 
Vdc = “37 V was used. Since the resulting DC electric field varies spatially, an accurate model 
of the trap must include the static potential z) which has the charge on electrode RF2 as 
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FIG. 2: Cross-section of the pseudopotential. The pseudopotential 'F (in units of eV) divided by the mass 
m (in units of kg) is plotted, so that the charge to mass ratio R found experimentally may be employed. 



FIG. 3: The pseudopotential (divided by mass) along the 2 ; direction at r = 0 (along the trap axis). A 
sixth-order fit is made to this potential. Writing the potential as /m{z — zq) = Yln=o ~ ^0)"'^ the fit 
parameters are cg = 9.87 x 10*^ eV/kg/mm®, C 5 = —1.11 x 10^ eV/kg/mm®, 04 = 5.37 x 10^ eV/kg/mm‘t, 
C 3 = 1.44x 10^ eV/kg/mm^, C 2 = 7.79x 10^ eV/kg/mm^, ci = —2.04 eV/kg/mm, and cq = —0.015 eV/kg. 
The range for this plot and fit is chosen to comfortably accommodate the range of motion of a trapped 
particle in our experiments. 
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its source. We use an effeetive foree model to describe the trap. We first define the effeetive 
pseudo-force as = — V^(r, z). In the following analysis, we eonsider only the z direction, 
as the ion is eonfined approximately to the z-axis. (We say “approximately” sinee, in our ex¬ 
periments, the amplitude of motion in the r direetion is on the order of one-tenth the amplitude 
in the z direction.) Considering this, the pseudo-foree is written as Frf = —d'^{z)/dz, where 
'1>(^) = ^(r = 0, z). The eorresponding foree from the applied DC eompensation voltage is then 
Fdc = —Qd^-Dc{z)/dz. The total foree is then given by 

F = Frf + -^DC ~ (5) 

where g is the aeceleration due to gravity. So that this expression may be written in terms of the 
eharge-to-mass ratio R of the ion, it is eonvenient to reeast it into an expression for the aeeeleration 
a{z) of the ion. In so doing, we obtain 

=-^^^\F{z)\^ FREbc- g, ( 6 ) 

where E{z) is the oseillating eleetrie field due to the potentials on RFl and RF3, and F^c is the 
static electric field due to the potential on RF2. In Fig.and its eaption we demonstrate the effeets 
of gravity and of the applied DC voltage on the aeeeleration experieneed by a trapped ion. Fig. 
shows the aeeeleration experieneed by a trapped ion near its equilibrium position, and a fit of that 
funetion is given. 


III. THE EXPERIMENT 

For this initial experiment with multipole traps, we use micron-sized granules of sugar whieh 
are held in air (inside an airtight eontainer). This is a suitable situation for this first set of measure¬ 
ments on SEMT’s sinee it is mueh simpler than a eomparable setup for atomie ions. However, it 
earries the disadvantage that the eharge and mass of the ion are not known in advanee. Therefore 
when we seek to mateh our experimental data with a model, we must adjust the model aeeording 
to various charge and mass values until the model eonverges to the experimental data. 

The trap is loaded by tapping eonfeetioner’s sugar from a plastie spatula onto the trap. At 
times, we use compressed air to blow the sugar off the trap surfaee. Normally, several partieles 
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FIG. 4: Plot of the acceleration experienced by a trapped ion as a function of position. The dotted (blue 
online) curve only accounts for the pseudo-force from the RF electrodes. The dash-dot (red online) curve 
shows the offset due to gravity; it is merely the blue curve translated down by 9.8 m/s^. The solid (green 
online) curve is the total acceleration, including the effect of the DC voltage of -37 V applied to RF2. The 
ion indeed is uncompensated (located below the point at which the RF electric field vanishes) when it is at 
the point at which the acceleration vanishes. However, when it is excited by an applied oscillating voltage, 
it is very nearly compensated at its maximum value of z. The total acceleration has a limiting value of 
-9.8 m/s^ as z —oo. 

are trapped and we eull the number to just one by manually adjusting the DC voltages. A single 
grain of sugar ean be stored for several hours; typieally, it is lost only when we shut the trap down. 
The ion seatters red laser light from a beam with a power of ~ 2 mW and a width of ~ 1 mm. The 
light is eolleeted using an imaging system with resolution 12 /rm and magniheation 3.4, whieh is 
ealibrated using a USAF 1951 resolution target. The seattered light is reeorded using a Thorlabs 
DCC3240C CMOS eamera. 

We measure the trajeetory of the ion as a funetion of time that results from an applied low- 
frequeney signal on one of the DC eleetrodes. At these low frequeneies, we simply use an op-amp 
adder eireuit to eombine the DC offset with the oseihating voltage, ealled the “tiekle.” This 
eombined signal is sent to an ampliher (Piezo Systems EPA-008-1) whieh has a gain of 20 in our 
frequeney range. The use of an op-amp eireuit has a niee side-beneht: the rails are hxed at +1 V 
and -7 V, whieh limits the input voltage to the ampliher; its maximum safe input is ±9 V. Tiekle 
frequeneies in the range of 1 to 50 Hz may be used, depending on the framerate of the eamera. The 
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FIG. 5: Plot of the total acceleration of a trapped ion, zoomed in on the true equilibrium position of the 
ion at Zmin = 1.4 mm. Although the potential has been altered by the need to compensate for the effect of 
gravity on the ion, the acceleration is still substantially nonlinear. A third-order polynomial fit is the solid 
line intersecting the computed data points. Writing the acceleration as a{z — Zmin) = Z]n=i {z — Zmin)^, 
the fit coefficients are = —87.7 m/s^/mm^, d 2 = 54.8 m/s^/mm^, and di = —13.4 m/s^/mm. 

maximum framerate depends on the size of the window, whieh itself depends on the amplitude 
of the ion’s motion in the trap. Typically, a window may be selected such that the framerate 
is approximately 60 which gives a Nyquist frequency of 30 Hz. Since a primary focus of 
the multipole trap measurements is frequency doubling, a nominal maximum tickle frequency is 
15 Hz. The data presented below were obtained using a measured tickle amplitude of 10.4 V. 

To obtain the data presented here, we record videos in avi format, typically consisting of 500 
individual frames. These videos are converted into a set of images using the software ffmpeg [l33l . 
which is called from Matlab, which is used to process the data. These images are then imported 
into Matlab, and the center of each ion is found. We use an algorithm which effectively calculates 
the “center of mass” for a given image, using the brightness of each pixel. These data in sequence 
form the position vs. time trajectory of the ion. Our method is very similar to that employed in 
the work of Ref. m, of which we became aware while our experiments were underway. 

Understanding the data requires comparing it with the model of the trap described above. This 
requires integrating the equations of motion for a single ion of a given R, then attempting to adjust 
R until a fit to the observed data is found. There is an additional term that depends on the mass 
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FIG. 6: Simulated (blue dashed line) and measured (red solid line) trajectory of an ion being excited at a 
frequency of 1.4 Hz by an oscillating voltage with amplitude 10.4 V applied to the electrode RF2. Our 
result values of i? = —2.54 mC/kg and /r/m = 674 s“^ are used in the simulation. The average height of 
the ion above the trap surface during the experiment is 1.4 mm (as shown in Fig. and Fig.[^. 

but not on the charge, which must be included in order to get an accurate simulation: the damping 
force from air. This force, denoted Fd, is described by 


Fd = -/iw, (7) 

where the constant /i quantifies the level of damping and v is the velocity of the particle in the 
direction. There are, in total, two parameters that must be fit, R and /x. If the ion were spherical 
in shape, we could use the result quoted in Ref. BUl, which states that 


/i = QtthdTp, ( 8 ) 

where /xd = 1-83 x 10“^ kg/m-s is the dynamic viscosity of air and is the radius of the particle. 
However, we lack the optical resolution to determine the shape of the particle; therefore we simply 
find the value of /x. 

In Fig. 1^ we present one example trajectory, for tickle frequency 1.4 Hz. The ion is clearly re¬ 
sponding to the sinusoidal tickle according to a nonsinusoidal waveform. The nonlinear properties 
of the trap may be most readily seen by Fourier transforming the position-vs-time data, revealing 
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FIG. 7: Fourier transform of the trajectory shown in Fig. The simulation is the blue dashed line, and the 
experimental data is the red solid line. 

the frequencies present in the trajectory. We find a strong peak at double the tickle frequency and 
a lesser peak at triple the frequency. An example of this Fourier transformed data is presented in 
Fig.|^ and Fig.j^is a two-dimensional plot of these data for several frequencies. 

For each trajectory, recorded at a different tickle frequency, we obtain a value for the charge- 
to-mass ratio R. The rms difference between the computed and measured trajectory is the error 
for each simulation, and the value of R which minimizes the error is found. A sample plot of this 
error as a function of simulated R is given in Fig.|^ For the particular ion used when the data were 
taken, we obtain a value of R = —2.54 it 0.01 mC/kg. This is a l-cx statistical error bar. There 
are, however, good reasons to believe that systematics dominate the error in the measurement, not 
the least of which is that important parameters such as the imaging system magnification and RF 
voltage are known to only two significant digits. We estimate the systematic error on R in the next 
section. Our measured value of R is well in line with previous measurements on charged dust 
particles OTI . The statistical error on the damping coefficient /i is greater because it affects the 
trajectory far less than R. For our result value of R, we obtain jj/m = 679 ± 46 kg/s. (Here, m 
is the particle mass.) This value is found in a similar way to the value for R, by minimizing the 
error between the experimental and simulated trajectories. The measurement of fi/m permits us 
to estimate the size of the trapped particles. They are almost certainly not spherical in shape, but if 
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FIG. 8: Two dimensional plot of the frequencies at which the ion responds versus the tickle frequency. The 
false-color plot is of the logarithm of the Fourier amplitude at each pair of frequencies. Frequency doubling 
is seen up to a frequency of about 6 Hz, after which it is difficult to see on this plot. One reason for this 
suppression of frequency doubling is the air damping, the force from which is velocity-dependent. 

they were, we eould invoke Eq. ([^ and use ps(/i/m)(47r/3)rp = where ps = 1590 kg/m^ 

is the density of suerose, to obtain a very reasonable value of 9 pm. This results in a partiele 
mass m of about 5 fg, and a eharge on the particle of Q = mR = 1 x 10“^^ C. These values 
are only estimates, but it appears in principle, for a particle whose drag is well-understood, to be 
possible to ascertain both the charge and the mass using our methods, not merely the ratio R. 


IV. DISCUSSION 

In this article, we have shown evidence that an ion trapped in our SEMT exhibits a nonlinear 
response to an applied signal. Although this has been seen in ion traps before, this experiment is 
unique in that the ion is close to the center of the trap during the experiment. The ion is actually 
close to being perfectly compensated when it is at its maximum height; we determine this by 
noting that the extent to which it is undergoing micromotion in the vertical direction is minimized 
at that point. Naturally, it would be ideal to repeat this experiment in vacuum, with an atomic ion 
of known charge and mass. One reason is that the effect of gravity on the atomic ion is negligible, 
and so a large compensation for the weight of the particle using an electric field sourced by RE2 
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FIG. 9: A typical plot of the error between the simulation and experiment as a function of simulated charge 
to mass ratio R. The error is calculated as the rms difference between the simulation and the experiment; 
therefore it has units of /rm. The red line connecting the data points is a fit of an eighth-order polynomial to 
the data. An eighth-order polynomial is used in our fitting routine so that wider ranges of R can be tested, 
but the value of R here does not differ substantially from that obtained by a quadratic fit. The minimum 
value of R is obtained by finding a local minimum of that fit function within the range of simulated R 
values used. 

would not be required, and we would therefore be able to see the true behavior of the multipole 
trap, without the quadratic terms introduced by RF2. Another aspect that would improve is that 
the damping of an atomic ion due to the laser cooling force would be easier to model than damping 
due to air, and we could therefore see a more direct confrontation between the experiment and the 
model. Nevertheless, with the experiment reported here, we see a clear qualitative signature of 
nonlinear dynamics, and find that our model yields a very reasonable value for the charge to mass 
ratio (i?) of the ion. We have also contrasted the observed multipole trap behavior with quadrupole 
trap behavior by trapping ions in quadrupole configuration (by applying the RF trapping voltage 
to RF2 while grounding RFl and RF3). We see no frequency upconversion in this case; the 
ion trajectory simply follows the drive. We also see the appearance of a secular frequency (the 
hallmark of a harmonic trap) in the Fourier transform of the trajectories. However, due to the 
presence of air, we do not see an increased amplitude of motion when the ion is driven at the 
secular frequency. Simulations verify that no resonant excitation should be observable, for our ion 
properties and other parameters. Future experiments in which we pump the system down should 
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allow us to observe resonant exeitation. As it stands, the shape of the potential itself depends 
strongly on the properties of the ion; therefore we cannot presently do a quantitative measurement 
of the size of each term in the potential. That said, there is no reason to doubt that the electrostatics 
used in the trap modeling are reliable. 

Since several of our experimental parameters are known to only two significant digits, we 
would expect systematic error on the measurement of R on the order of a few percent. We estimate 
the systematic error in the following way. We note that the shape of the trajectory alone should be 
sufficient to calculate R, without regard to the absolute amplitude of the trajectory. This is because 
the trap acceleration function is both non-linear and non-periodic; there is only one position along 
the curve that can generate a particular response trajectory from the ion. Therefore, we repeat 
the fitting process, but normalizing the trajectories each time (to have a maximum value 0.5 and 
minimum value -0.5). For each dataset, the value of R that minimizes the error is found. Using 
this approach, we obtain R = —2.69 ± 0.04 mC/kg. This is an error of 6% compared with the 
other measurement, in line with our expectations. We have confined a few other sugar particles 
in our trap, and done a rough estimate of their R values using the methods outlined in this paper. 
The mean value is around 7 mC/kg with a standard deviation (over six measurements) of 5 mC/kg. 
This ought not be surprising, since they are charged triboelectrically. 

Since these measurements represent a new way (to our knowledge) of determining R for some 
substance, it is natural to ask whether they may be adapted into a new type of mass spectrome¬ 
try. Perhaps, since only the shape of the trajectory needs to be known (rather than the absolute 
amplitude), the excitation of the motion of the unknown ion could be detected electronically. Fur¬ 
thermore, SEMT’s have several other desirable properties. SEMT’s, in contrast to quadrupole 
traps, appear to be able to simultaneously confine particles of widely different R values. In addi¬ 
tion, the depth of the trap, as commonly calculated for quadrupole traps, is significantly greater 
than that of quadrupole surface traps of similar scale [|24l . This leads us to the idea that an array 
of SEMT’s could work as a mass analyzer, providing a broad spectrum spanning many R values 
simultaneously. The traps would need to be very small, and multiple traps would need to be op¬ 
erated in parallel, but that is established technology; a microfabricated array of quadrupole ion 
traps for mass spectrometry has already been demonstrated ll3^ . Eurther work on optimization 
of the trap geometry, as well as mapping of the stability diagram of the traps, will help to resolve 
whether this is a possibility. Successful operation would require marked improvement in system- 
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atic sources of error, whieh would likely be possible when trapping atomie or moleeular ions, 
espeeially if the trap ean be ealibrated using a partiele of known R. 
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